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ABSTRACT 

Recently, the detection of discrete features in the X-ray afterglow spectra of 
GRB970508 and GRB970828 was reported. The most natural interpretation of these 
features is that they are redshifted Fe K emission complexes. The identification of 
the line emission mechanism has drastic implications for the inferred mass of radiat- 
ing material, and hence the nature of the burst site. X-ray spectroscopy provides a 
direct observational constraint on these properties of gamma-ray bursters. We briefly 
discuss how these constraints arise, in the context of an application to the spectrum of 
GRB970508. 

Subject headings: gamma rays: bursts — X-rays: general — techniques: spectroscopic 
— line: formation 



1. Introduction 



The detection of a discrete spectral feature in the X-ray afterglow spectrum of GRB970508 was 
recently reported by Piro et al. ( 1999a, b). A similar feature in the afterglow of GRB970828 was 
reported by Yoshida et al. (1999). The redshift of the host galaxy of GRB970508 was determined to 
be z = 0.835 (Metzger et al. 1997; Bloom et al. 1998). The apparent energy of the discrete feature 
in the X-ray spectrum is consistent with the redshifted energies of the n = 2 — 1 transitions in all 
possible charge states of Fe, i.e., from Fe Ka at 6.4 keV, to H-like Fe Lya at 6.95 keV. Therefore, 
none of the possible line excitation mechanisms (fluorescence, collisional excitation in hot gas, or 
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recombination in photoionized gas or a transient collisional plasma) is currently ruled out on the 
basis of the measured line energy alone. 

The conditions required for each of these emission mechanisms (i.e. the presence of a given 
charge state of Fe) are all still compatible with the constraints on ionization parameter and gas 
temperature that can be derived from the fact that the line source must be located no more than 
about 1 light day from the ionizing source, and from the shape of the ionizing spectrum. All three 
mechanisms require high gas density (n e ^ 10 11 cm~ 3 ) for bound Fe to exist at all this close to the 
burst site. The presence of a large mass of dense gas close to the burst has led to the suggestion that 
a merging neutron star- neutron star binary is an unlikely site for this burst, and that a scenario 
involving some sort of stellar collapse is more likely (Piro et al. 1999b). In that case, the line source 
may be associated with debris from the stellar collapse. 

Resolved X-ray spectroscopy can of course distinguish between the various line emission mech- 
anisms, and lead to the correct characterization of the conditions in the source and the implied 
mass of radiating Fe. In the following, we will briefly recapitulate the constraints imposed by 
proximity and the energy budget, and then draw attention to an example of direct spectroscopic 
analysis of the X-ray spectrum: the fact that the measured centroid position, and to some extent 
even the shape, of the X-ray emission feature in the afterglow of GRB970508 already imply that 
recombination in afterglow-photoionized gas is unlikely to be the emission mechanism for this burst. 

Our analysis applies to an optically thin, homogeneous medium. Complications arising from 
radiative transfer effects or external heat input by a relativistic shock may conspire to alter the 
shape of the spectrum, perhaps to the point that no definitive conclusions can be drawn from an 
intantaneous, low sensitivity spectrum. Such effects have been calculated by Weth et al. (1999) 
for the case of photoionization in the context of two specific geometries for the burst site, and by 
Vietri (1999) for a scenario in which the relativistic shock heats the medium shortly after the onset 
of the afterglow. With better data, it may hopefully be possible to reverse the argument, and infer 
the geometry, thermal history, and abundances from spectroscopy. 

2. Ionization and Thermal Conditions in the GRB970508 X-ray Line Source 

The discrete feature appears in the BeppoSAX Medium Energy Concentrator Spectrometer 
(MECS) spectrum immediately before a sudden rise in the afterglow flux, beginning at ~ 6 x 10 4 
sec after the burst, and disappears thereafter (Piro et al. 1999b). If one attributes the feature to 
emission by Fe, one requires bound Fe to be present close to the source, which emits a very large 
flux of ionizing photons in the afterglow (L\ R y( j-io keV ~ 10 51 (t/l sec) -11 erg s _1 ; here, t is the 
time since the burst [Piro et al. 1999b]). 

The prompt burst and afterglow are intense enough that the photoionization timescale, £i on , is 
extremely short compared to the burst duration and the time during which the emission feature is 
visible. Specifically, the (inverse of the) ionization timescale for Fe in the afterglow radiation field 
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I = j dE F(E) • E~ l a{E) = (1) 
1.4 x 10 4 (t/1 sec)" 1 ' 1 rfg 2 sec" 1 , 



with F(E) the ionizing flux at the line source, a(E) « 6 x 1CT 18 Z~ 2 (x/E) 3 cm 2 the K-shell 
photoionization cross section for neutral atoms of nuclear charge Z, \ the ionization potential, 
and ri6 the distance to the burst site in units 10 16 cm. We have assumed a E~ 2 photon number 
spectrum to calculate F{E) from the 2 — 10 keV afterglow flux quoted by Piro et al. (1998), 
converted to luminosity assuming Hq = 75 km s _1 Mpc _1 and qo = 1/2, and we extrapolate the 
luminosity to 1 MeV. Ionization is essentially instantaneous. 

The recombination timescale is given by i rcc ~ (n e a(T e )) _1 , where a(T e ) is the recombination 
coefficient, n e the electron density, and T e the electron temperature. Since a depends on the 
temperature, we first need to consider the thermal evolution of the source subject to irradiation by 
the prompt burst and afterglow radiation fields. 

Given sufficient interaction time, the gas will relax to the Compton temperature, Tq, at which 
the Compton heating and cooling rates equal each other, and which is therefore determined only 
by the shape of the ionizing spectrum. In the non-relativistic limit (both the photon energies E 
and the electron energies kT e <C m e c 2 ), the Compton temperature is given by 

1 f dE E F(E) 

kTc = 4 Si EF(E) (2) 

(Ross 1979). As either the photon energies or the gas temperature approach the electron rest energy, 
the correct relativistic scattering cross section should be used, and the energy exchange between 
photons and electrons should be calculated to higher order than linear in E/m e c 2 and kT e /m e c 2 . 
In order to obtain a very rough estimate, we will assume here that we can mimic the effects of these 
modifications by simply cutting off the integrals at photon energies of order 1 MeV. If we assume a 
prompt burst spectrum of the form F(E) cx E 1 ! 2 at low energies, and F(E) cx E~ 3 / 2 above a break 
energy of order 1 MeV (Conners et al. 1998), the Compton temperature becomes approximately 
kT c (MeV) ~ (3/16)(£ max /l MeV) 1 / 2 , with E max the cutoff energy. For E max in the few MeV 
range, the Compton temperature is therefore in the several hundred keV range. The timescale for 
the Compton interactions to equilibrate should be of order (Rybicki Sz Lightman 1979) 

tcompton = , ~ 7 x 10 n 10 T g sec (3) 

with <tt the Thomson cross section. Unless the density is much higher than 10 10 cm~ 3 , this timescale 
is much longer than the burst timescale, which implies that the Compton temperature will not be 
reached. Moreover, the average photon energy in the burst spectrum is of order 1 MeV, which 
produces mildly relativistic electrons. The stopping timescale of these fast particles on stationary 
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electrons is (e.g. Longair 1992) 



t, 




= 1.3 x 10 3 nr r 
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( 7 /2) ((ln^/lO)" 1 sec, 
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(with 7 the electron Lorentz factor, and In A the Coulomb logarithm) which leads to the interesting 
conclusion that the plasma may not be in equilibrium at the end of the burst. 

In any case, as the (much softer) afterglow begins, the plasma will rapidly cool by inverse 
Compton scattering, bremsstrahlung, and, at low temperatures, atomic emission. For an initial 
gas temperature of order 10 8 K or hotter, the Comptonization timescale is of the order of, or 
shorter than the duration of the afterglow, and equilibrium is established at the afterglow Compton 
temperature, which for an E~ 2 photon spectrum is equal to kTc = (l/4)E max / hi(.E' max /I£ m i n ), with 
.Emit! the lowest photon energy in the afterglow spectrum; E max is the upper end of the integration 
range, -Emax ~ 1 MeV. Assuming the afterglow spectrum flattens below E m \ D ^ 1 Ryd (Wijers & 
Galama 1999), we find Tq < 2 y. 10 8 K. Kallman and McCray (1979) in their numerical models 
for X-ray photoionized nebulae, actually find that the electron temperature saturates at a lower 
temperature of approximately Tq ~ 10 7 K, due to bremsstrahlung cooling, so most likely the gas 
is cooler than 10 8 K in the afterglow. 

The recombination coefficient as a function of temperature for recombination onto a bare 
nucleus is given by Seaton (1959); a powerlaw fit to the coefficient is a ~ 1.26 x 10~ 6 T~ ' 75 cm 3 
s _1 (for Z = 26), accurate to 20% in the range T e = 10 5 — 2 x 10 s K. Thus, the recombination 
timescale is 



We conclude that at sufficiently high densities (nn <; 1), the line source will be relatively cool, and 
close to ionization equilibrium at all times during the afterglow. 

Assuming ionization equilibrium, we can derive a constraint on the density from the require- 
ment that bound Fe be present in the line source. For Fe to recombine to the H-like charge state 
requires an ionization parameter in equilibrium of £ = L/nr 2 ^ 10 4 , with L the ionizing luminos- 
ity, n the particle density, and r the distance to the ionizing source. This value of the ionization 
parameter applies to an E~ 2 ionizing photon spectrum (Kallman &i McCray 1982, their model 7). 
Inserting numbers and extrapolating the ionizing luminosity to 1 MeV, one infers a lower limit to 
the density at a distance of 6 x 10 4 It sec from the source, of 



Similarly, for Fe to be recombined to Li-like or less ionized, necessary for fluorescence to operate, 
£ 1000 in equilibrium, yielding a ten times higher density than the above estimate. 



t rec = 1.4 n n 4 T 7 - 75 sec. 



(5) 



n Z 1.7 x 10 11 (t/6 x 10 4 sec)- L1 (r/6 x 10 4 It sec)" 2 cm" 3 



(6) 
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3. X-ray Spectroscopy 

From the previous we conclude that if the interaction with the afterglow radiation is the only 
source of heating in the line emitting gas, then, later in the afterglow, the source is likely to be 
dense 10 11 cm~ 3 ) and relatively cool. Line emission can be excited both by fluorescence and 
cascading following recombination. Only if there is an additional source of heat, sufficient to heat 
the gas to T ^ 10 8 K, will collisional excitation play a role in driving Fe n = 2 — 1 lines. 

Obviously, X-ray spectroscopy of the discrete emission strongly distinguishes between these 
various possible emission mechanisms. If one is confident of the redshift of the emission complex, 
fluorescence can be distinguished from emission by highly ionized Fe, simply from the apparent 
energy of the emission. Unfortunately, the resolution of the MECS, and the statistical quality of 
the spectra under discussion are not sufficient to allow this distinction to be made. 

Recombination and collisional excitation, however, can readily be distinguished in the follow- 
ing way. For a given degree of ionization, a photoionized medium in equilibrium is much cooler 
than the corresponding collisional plasma, in which the electron temperature has to be of order 
the ionization potential in order to support a given charge state. The low electron temperature in 
photoionization equilibrium has a dramatic spectroscopic signature. Since most of the free electrons 
have small kinetic energies compared to the ionization potential, photons resulting from radiative 
recombination will have a narrow energy distribution, bunched up just above the ionization po- 
tential, with a typical width of order AE ~ kT e «x. In a collisional plasma on the other hand 
(kT e ~ %), the recombination rate is much reduced, and in addition, the recombination photons 
are spread out over a wide energy range above the ionization potential. The narrow radiative re- 
combination continuum (RRC) in a photoionized source is easy to detect, because it contains an 
integrated photon flux comparable to that in the n = 2 — 1 discrete transition; in fact, for H-like 
Fe, the ratio between the fluxes in the RRC to that in the Lya line is approximately equal to 
0.93 (kT e /l keV) - 21 . So far, narrow RRC's have only been seen in the spectrum of the massive 
binary Cyg X-3 (Liedahl & Paerels 1996; Kawashima & Kitamoto 1996). 



4. Application to the Spectrum of GRB980508 

We use the dataset shown by Piro et al. (1999b; their dataset la). At the resolution of the 
MECS (AE ~ 340 eV at 3 keV), the RRC is barely resolvable if the electron temperature is of order 
or less than 1 keV. For simplicity we use a model for the emission spectrum of cool, photoionized 
gas consisting of a narrow line plus exponentially decaying RRC, to represent the Lya emission line 
and RRC at 6.95 and 9.28 keV, respectively, with equal photon fluxes. The He-like recombination 
spectrum would look much the same at the MECS resolution. To this, we added a simple power 
law with absorption by neutral gas, to represent the continuum. We let the redshift float, as well 
as the continuum parameters. 
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By fixing the electron temperature (i.e., the width of the RRC) and fitting for the remaining 
free parameters we obtain the minimum x 2 values displayed in Figure 1. At either very low T e 
(kT e ^ 0.2 keV) or very high T e (kT e <; 4 keV) we obtain x 2 values of order Xmin ~ ^ f° r ^ 
degrees of freedom. In themselves these values are acceptable (given the small number of degrees 
of freedom). However, at low temperatures the model has a best fitting redshift of z > 1.6 (for 
Hydrogenic Fe; 68% confidence for one parameter of interest only). At high temperatures, the 
RRC becomes essentially undetectable and the model consists effectively of a single narrow line. 
This model will of course fit the data, but at these temperatures the plasma would be closer to 
collisional equilibrium, which is not consistent with the assumption we set out to test. Finally, at 
kT e ~ 1 keV (which is about the temperature for a source in photoionization equilibrium with the 
afterglow), we find Xmin ~ 12, which indicates a poor fit (Xmin about two standard deviations away 
from the expected Xmin = 5)- By itself, this is probably not enough to confidently reject the model, 
and the fit can be improved by changing the temperature somewhat. But the implied redshift is 
still z ~ 1.20, higher than for a single line because we have an additional emission component at 
4/3 times the line energy. Unless we are willing to entertain the possibility that the line source 
is actually behind the galaxy at z = 0.835, we conclude that the emission spectrum from a cool 
photoionized source is incompatible with the measured spectrum. In Figure 1, we have indicated 
the (68% confidence) limits on the redshift as a function of the assumed temperature of the gas, 
assuming H-like emission. At high temperature, the contrast in the RRC becomes very small, and 
the spectrum is dominated by a single line. In this limiting case, the implied redshift approaches 
~ 1.05 ±0.05 (68% confidence), appropriate for H-like Fe Lya. Widening the confidence interval to 
90% confidence will produce agreement with the optical redshift of z = 0.835 (Piro et al. 1999b). 

Figure 2 displays the data and the lowest— x 2 recombination model with kT e = 1 keV; the values 
for the other spectral parameters are: power law photon index —3.1, column density Nh = 7.1 x 10 21 
cm~ 2 . This continuum spectral shape appears to be quite a bit steeper than the best fit given by 
Piro k, al. 1999b, but in fact, the uncertainty on the index is large (w 1.1 for one parameter of 
interest, at la), and so the continuum shapes are in fact consistent. 

5. Conclusion 

We find that Fe emission from a cool photoionized source in equilibrium with the afterglow 
cannot account for the X-ray spectrum of GRB970508. Fitting a spectral model for such a source 
either implies a large redshift, or indicates such high electron temperatures that the source would 
be closer to collisional ionization equilibrium. 

This would imply that the discrete feature in the spectrum of GRB970508 arises from either 
Fe fluorescence or from collisional excitation, each of which implies a very different estimate for 
the total mass of radiating Fe (M Fe ~ 2.6 x 10" 4 M© for fluorescence, M Fe ~ 7 x 10~ 2 M© for 
collisional excitation, which corresponds to M to tai ~ 8 x 1O _2 M (^IpeMFe,©) -1 and M to tai ~ 
28M Q (ylpe/^-Fe,©) -1 , respectively (see, for instance, Meszaros & Rees [1998], Piro et al. [1999b], 
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Lazzati et al. [1999]); here, Ap e is the abundance of Fe). 

With higher spectral resolution and sensitivity, we may be able to distinguish spectroscopically 
between a fluorescent and a collisional spectrum (a recombination spectrum is readily identified 
from the presence of the RRC). Working from the K-shell spectrum of a single ionization stage of 
Fe alone, one might use the n = 2 — 1 and n = 3 — 1 transitions, which are resolved at E/AE ~ 
10. Intensity ratios do not effectively discriminate between collisional excitation and fluorescence: 
/(Ly/3)/I(Lya) = 0.10 — 0.15 (Hydrogenic Fe) for electron temperatures between 3 x 10 7 and 
10 9 K (Mewe, Gronenschild, & van den Oord 1985), while for neutral Fe I(K(3)/I(Ka) = 0.13 
(Kaastra & Mewe 1993). Instead, with somewhat higher resolution, one might use the fact that 
the ratio of the transition energies depends on the charge state: E(Ly(3)/E(Lya) = 1.18, while 
E(Kf3)/E(Ka) = 1.10. 

At higher resolution, one would resolve the ls — 2pi/2 an d Is — 2p 3 / 2 transitions (E/AE > 500). 
The ratios of the energies are E(Lya 2 )/E(Lya 1 ) = 1.0030 and E(Ka 2 )/E(Ka 1 ) = 1.0020, and a 
resolving power of 1000 is required to distinguish between these two. 

A collisional source at kT < 15 keV will also show emission from He-like Fe, and the simul- 
taneous appearance of both the H- and the He-like charge states readily indicates a high degree 
of ionization of the source. For temperatures in the range 3 < kT < 15 keV, the He-like emission 
dominates, which could still be identified as such if the characteristic "triplet" structure of the 
n = 2 — 1 lines is resolved (E/AE > 300). The highest resolution available at energies above 2 keV 
is provided by the High Energy Transmission Grating Spectrometer on Chandra (Canizares et al. 
2000). The Fe K spectrum will be resolved at E/AE s=s 150 (1 + z), which would at least enable 
charge-state spectroscopy at moderately large redshift. 

Finally, given sufficient sensitivity at low energies (and sufficiently low absorption in the 
source), we may of course detect emission from lower- Z elements, either collisional or fluorescent, 
which would greatly facilitate the spectroscopic diagnosis of the excitation mechanism. With the 
grating spectrometers on Chandra and XMM, standard detailed plasma diagnostics can be invoked, 
at resolving powers of E/AE ~ 100 — 1500 at E ^ 2 keV, depending on photon energy (Canizares 
et al. 2000, Brinkman et al. 2000, Brinkman et al. 1998). Clearly, it would be very interesting to 
pursue more sensitive X-ray spectroscopy of gamma ray burst afterglows. 

We gratefully acknowledge discussions with Marten van Kerkwijk. D. A. L. was supported in 
part by a NASA Long Term Space Astrophysics Program grant (LTSA S-92654-F). Work at LLNL 
was performed under the auspices of the U. S. Department of Energy, Contract No. W-7405-Eng-48. 
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Figure Captions: 

Fig. 1. — Minimum \ 2 and best fitting redshift for fitting a spectral model appropriate for emission 
from afterglow-photoionized gas to the BeppoSAX MECS+LECS spectrum of GRB970508, as a 
function of electron temperature kT e . The grey band represents the la confidence volume for the 
X-ray redshift; the dotted lines indicate the la contours for the best-fitting X-ray redshift assuming 
H-like Lya line emission only (no RRC). 

Fig. 2. — Best fitting photoionization equilibrium emission model, for an assumed electron temper- 
ature of kT e = 1 keV. Solid triangular datapoints are LECS fluxes, solid circles are MECS fluxes. 
This model implies a source redshift of z = 1.20 for H-like Fe emission. Lower panel displays the 
post-fit residuals. 
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